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Recent developments on studies of transport through quantum dots obtained by applying the time-dependent 
density matrix renormalization group method are summarized. Some new aspects of Kondo physics which 
appear in nonequilibrium steady states are discussed both for the single dot case and for the serially coupled 
double-quantum-dot case. 
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1 Introduction 

Advancements in mi crof abdication technology provide us opportunities to investigate transport properties 
through quantum dots. Stimulated by possible application to quantum computation, number of interesting 
systems are now being produced by arranging multiple quantum dots in various geometries. 

In the simplest case, namely a single quantum dot system, it was predicted [[I] and experimentally 
observed j2j that the linear conductance is enhanced at low temperatures up to the unitary limit by the 
Kondo effect when odd number of electrons are accommodated in the dot. 

To measure a current through quantum dots a finite bias voltage is applied between source and drain 
electrodes. Therefore, in principle, we always confront the problem of nonlinear transport. Away from the 
linear response regime, the differential conductance for single dot system is quickly suppressed reflecting 
the sharp Kondo resonance peak at the Fermi energy. The nonlinear transport through quantum dots defines 
a prototypical problem for electron correlation under a nonequilibrium condition. Various theoretical ap- 
proaches have been applied to this problem but our understanding of the problem is still far from complete. 
Therefore it is desirable to obtain unbiased results by using some computational approaches. 

Recently there has been progress in computational physics for strongly correlated electron systems from 
two complementary approaches. One is from the infinite dimension initiated by Metzner and Vollhardt 
which was later developed into the dynamical mean field theory (DMFT) (4]. It is interesting to note that 
the Kondo effect is relevant to the DMFT since in the infinite dimension the self-energy becomes local, 
which is determined by similar physics as the Kondo effect. The other is the density matrix renormalization 
group method developed by White |5 |, which is an approach from one dimension. Among many numerical 
techniques which have been extended to nonequilibrium situations, the time-dependent density matrix 
renormalization group method (TdDMRG) [[6] has proven to be a powerful tool to tackle the nonlinear 
transport in strongly correlated electron systems, such as the single quantum dot system JH] [9] [TT] [TOj and 
the interacting resonant level model |[T2| [T3l . In this paper we discuss our recent activities to understand 
the nonequilibrium transport through quantum dots by using the TdDMRG. We will briefly summarize our 
results concerning the single quantum dot J5] [TT] and then show recent results on the serial double-dot 
system. 
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2 Single Quantum Dot System 

We consider a system consisting of a quantum dot and two leads. We focus on the situation where the level 
spacing in the quantum dot is large compared to the other energy scales in the system, so that only one 
level in the dot is active and the others are effectively frozen. Concerning the leads we assume that single 
channel is relevant to transport, which can be described by the 1-D tight-binding model. The Hamiltonian 
for the system is expressed as follows: 

H ( T ) = - t H^i+la + h.C.) -i^^(4c i+1(T + h.C.) 

i< — l <x z>0 cr 

-t'Yl [( C -l- C 0- + k - C -) + ( C CT Cla + h.C.)\ 

~ | E c o« c °° + ^ot^o+cot + \ 6 ( T ) {Nl - N R ) , (1) 

a 

where r is the time variable, the quantum dot is located at the Oth site and c\ a creates an electron with spin 
a at ith site. t,t' ,U and V are the hopping amplitude in leads, the one between the dot and the neighboring 
sites, the Coulomb repulsion energy and the voltage drop between the two leads, respectively. Here we 
assume the electron-hole symmetry, which means that the one-particle energy at the dot site is fixed as 
— U/2, and the inversion symmetry at the dot site except for the voltage term. Nl is the sum of the number 
operators in the left lead and Nr in the right lead. In the following we focus on the half-filled case. 

In the Keldysh formalism, one starts with an equilibrium state of an unperturbed Hamiltonian. One then 
turns on the perturbation term adiabatically, and gets a nonequilibrium steady state after a relaxation time. 
Although there are several choices for dividing the Hamiltonian into unperturbed and perturbed terms, the 
nonequilibrium steady states for the full Hamiltonian are expected to be independent of the choice. For 
the present problem we take a time-dependent voltage term. 9(r) is a smoothed step function which is 
introduced to simulate adiabatic switching-on of the voltage. Throughout this paper we set 6{t) as 




The reason for our choice of the time dependence of the Hamiltonian is that, if t' is chosen as the per- 
turbation, there is only one nonzero eigenvalue of the reduced density matrix of the lead parts, i.e. the 
entanglement entropy between dot and leads are zero for the initial states. Then the usual DMRG scheme 
to obtain the optimal basis set, which will be necessary for the subsequent TdDMRG calculation, cannot 
be used directly. 

By the TdDMRG method we can compute the time evolution of the initial wave function. In order to 
obtain the current at each time we can simply take the expectation values as follows: 

J(t) = -<V(r)|eJVi|V(r)) = — ^(V>(r)|(4 CT C-i CT - fc.c)|^(r)>. (3) 

cr 

Note that one can define the current operators for every bond in the system. When the system reaches a 
steady state, the currents stay constant at a certain value because of the definition of the steady states and 
the charge conservation. 

The steady current in the thermodynamic limit can be calculated based on the Keldysh formalism as 

OH 

J(V) = £ E ^f^fe B) Im ' (4) 
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Fig. 1 (Left) Current-voltage characteristics of the Anderson model Eq.^TJ. (Right) Differential conductance obtained 
by interpolating J(V). We omit the data where numerical differentiation is so delicate that it can produce an unphysical 
peak structure in G(V). Adopted from Ref. llll . 



where Tl r{cj) is the resonance width due to the hybridization to the left and right leads which will be 
assumed to be the same in this paper, and g^ ot (uj) the retarded Green's function at the dot site. When 
{7 = 0, the problem is simplified to the one -particle level. In this case Eq.© can be exactly evaluated. 
Alternatively, one can numerically diagonalize the one-particle Hamiltonian and explicitly calculate the 
time-evolution operator for a finite size system. Errors included in TdDMRG results can be estimated by 
comparing with results obtained from these methods. 

Since our system is finite we cannot reach to real steady states. However, it is possible to evaluate steady 
currents from the quasi-steady states which are defined in the time domain after the initial relaxation time 
and before the arrival of the wave front reflected at the boundary of the system. 

Quasi-steady states can be realized by the TdDMRG calculation for U/t < 5. The fixed value i'/i = 0.6 
and the range for U correspond to the Wilson ratio 1 < R\y < 1.97. Then we can extract the nonlinear 
current-voltage characteristics from the quasi-steady regions, which are summarized in FigsfTI IfTTI . In 
the noninteracting case the TdDMRG results perfectly agree with the exact analytic results for L — > oc. 
Moreover, the linear conductance for every U corresponds to the prediction of the Fermi liquid theory JT] , 



dJ_ 
W 



2e 2 

= -T- (5) 

v=o n 



In the differential conductance G(V) we can see a substantial difference in the width of the zero bias 
peak depending on U which is a reflection of the sharp Kondo peak in the local density of states at the dot 
site. Therefore we can say that the TdDMRG calculation catches the essence of the Kondo effect of the 
quantum dot system out of equilibrium in the parameter space investigated. 



3 Serial Double-Quantum-Dot System 

Now we proceed to the serially coupled double-quantum-dot system |[T5l[T6l[r7l[T8l , each connected with 
a lead, and discuss transport properties in the system driven by a DC voltage. From a theoretical point of 
view the system can be modeled by the two impurity Anderson model. 

The model is depicted in Fig|2]and it is straightforward to write the Hamiltonian explicitly in the similar 
form as Eq.dTJ. The only new additional parameter is the interdot hopping, —t". In what follows we 
discuss the current between the two dots, 

jpf'i 

J M = xE^WI(^ - h.c.)Mr)), (6) 

a 
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where ct (c£ CT ) creates an electron at the left (right) dot. 

There are several theoretical studies on the double-quantum-dot system. For the linear transport, the 
slave boson mean field theory (SBMFT) [19 ,'20'j and the numerical renormalization group (NRG) method 
ETl have been applied to the problem and predicted that the linear conductance has a sharp peak which 
reaches the perfect conductance 2e 2 /h as a function of the ratio of the interdot hopping amplitude t" to 
the resonance width at the Fermi energy V, which is given by T = t' 2 /t for the present model. The peak 
is located at J e g ~ X*., where J c s = denotes the antiferromagnetic exchange coupling between the 
two dots, and the Kondo temperature at t" = 0. 

This behavior can be understood in connection with formation and competition of the various types of 
singlets in the system as explained in Ref. ED . The conductance formula is expressed by 11221 

G=^sm 2 (6 e -5 ), (7) 

where S e and S a are the scattering phase shifts for even and odd channels, respectively. The even and odd 
orbitals of the two dots are defined as c\ a = -j= (c\ a + cJ CT J and c oa = ^cj a — cJ^J . When each 
impurity site contains one electron, S e + S a = it follows from the Friedel sum rule. Then, the following 
three characteristic states of the two impurity Anderson model appear depending on t", T and U: 

(i) For small t" /V such that J c ft <C TJJ-, the Kondo singlet state is formed on each dot with its adjacent 
lead. The phase shifts for this state are S e <~ S a ~ n/2. 

(ii) In the intermediate t" /T region with J e g 3> T^, the local spins on both dots are coupled as a spin 
singlet state, ^(cj^c^ - c^c^)|0): apair spin singlet. S e ~ it, 6 ~ 0. 

(iii) When t" > C//4 the doubly occupied state of bonding orbital of the two dots, cL-C^O), becomes the 
most stable: a bonding singlet. S e ~ n, S a ~ 0. 

The above three states are schematically shown in Figf3]and their relationship is displayed in Figf2](b). 
For the ground state (i) the two dots are effectively decoupled and similarly for (ii) and (iii) the two dots 
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Fig. 2 (Left) Model of the serial double-quantum-dot system. (Right) Schematic diagram for the characteristic ground 
states (i), (ii) and (iii) on the t" — U plane. The curves represent crossover between the states. 




Fig. 3 (Left) Schematic figure of the two-Kondo-singlets state. (Middle) The pair spin singlet state. (Right) The 
bonding singlet state. 
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Fig. 4 Current-voltage characteristics of the serial double-quantum-dot system for (7 = obtained by using exact 
diagonalization. Results for t' /t — 0.6 and t" /t — 0.6 are obtained for system of L = 81 sites and the others 
are for L = 121. Arrows indicate the positions of the inflection points, that is, the peak positions in the differential 
conductance d.J{V)/dV . Curves without associated arrows have inflection points at V = 0. 



can be viewed as detached from the two leads, resulting in small conductance for all of the above states. 
A sharp peak with height 2e 2 / h is found in the crossover region between (i) and (ii), namely J c ff « T^-. 
The continuous changes of 5 e and S from (i) to (ii) (or directly to (iii) for relatively small U) with keeping 
the Freedel sum rule ensure that Eq.(|7) has the maximum which reaches the perfect conductance. Just at 
the peak the system forms a coherent superposition of the two-Kondo singlets on the dots. As U increases 
or r decreases, in other words the electron correlation becomes strong, Tj ( rapidly decreases and the peak 
position moves to smaller values of t" . In addition the width of the peak becomes narrower, representing 
the decrease of the Kondo temperature. These intriguing behaviors can be attributed to the quasi-quantum 
criticality of the two impurity Kondo problem. 

Let us first discuss the current-voltage characteristics in the noninteracting case. Since the parameters 
for the two dots are identical and the energy levels are set to be particle-hole symmetric, each dot contains 
exactly one electron. Then for U = the linear conductance is given by 

h [i + (f /r) 2 ] 2 

G becomes the perfect conductance 2e 2 /h at the point t" /V = 1. From FigJU the linear conductance for 
(t'/t,t"/t) = (0.6,0.4) and (0.8,0.6) is almost the perfect one (1.98e 2 //i and 1.99e 2 //i), since t"/T = 
1.11, 0.94 is close to unity. 

The LDOS for the even and odd states have peaks at ±<" with widths ~ T. The electrons moving 
through the two quantum dots, which are inside the bias window — eV/2 < ui < eV/2, use these peaks in 
the resonant tunneling. When T is large, the two peaks merge into a single peak located at oj = 0. On the 
other hand for t" > T, the differential conductance G(V) = dJ ^y has its maximum when the peaks in 
the LDOS come into the bias window. Concerning the results shown in Fig0J for t' jt = 0.6 and 0.8 the 
maximum of the slope of J(V) is located at V = 0, reflecting the single peak structure of the LDOS. On 
the other hand for t'/t = 0.4 we observe the split peaks of G(V) at eV = ±2t" . 

Now we consider the effects of the Coulomb interaction on the currents. First we study the linear 
conductance and it is summarized in Fig|5] For each U the conductance as a function of t" /V forms a peak 
near t"/T = 1 and their heights reach almost 2e 2 /h JH1 |T9l [20l EU . While on the right side of the peak the 
conductance decreases as U increases, it increases on the left side, reflecting the shift of the peak toward 
smaller t"/T, in accordance with the NRG results ifZTl . 

Several typical examples of the nonlinear currents of the double-quantum-dot system are shown in 
Fig0 [23], together with the differential conductance G(V) = and the expectation values of the 
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Fig. 5 Linear conductance of the serially coupled double-quantum-dot system obtained from J(V) /V\ e v/t=o.i- 
The sets of data are for (t'/t, t"/t) = (1, 0.4), (1, 0.6), (0.8, 0.4), (0.8, 0.6), (0.6, 0.4), (0.6, 0.6), (0.4, 0.4) and 
(0.4, 0.6) from left to right. 



spin correlation operator (Si ■ S r ). Note that the results for G(V) have larger errors than the others since 
they are obtained via interpolation and differentiation of the data points of J{V). 

The results for (t'/t, t" jt) = (0.8, 0.6), not shown in this paper H231 . are close to the linear conductance 
peak at t" /T ~ 1. The linear conductance is thus near 2e 2 /h and the system as a whole forms a coherent 
singlet state. Out of the linear response regime we find a substantial decrease of the current due to the 
Coulomb interaction. This can be attributed to the coherent singlet state whose weight of the LDOS builds 
up in a narrow region around the Fermi energy, similarly to the case of the single dot. 

Next we discuss the results on the right side of the linear conductance peak at t"/T ~ 1, that is, the 
results for (t'/t, t"/t) ~ (0.4, 0.6), left panels in Fig|6l In this case the ground state has large amplitudes 
on the states (ii) and (iii), which can be confirmed by relatively large values of \(Si ■ S r )\. The splitting of 
the zero bias peak can be seen in G(V ) in this case. The position of the peak eVp ea k moves to smaller V as 
increasing U and corresponds to the position where \(Si ■ S r )\ varies significantly. This can be explained 
as follows: as V becomes comparable to V^, e ak the local electronic state is modified in the steady state. 
Hence it is confirmed that the many -body bonding and antibonding states are located at ±eVp e ak/2 in 
the LDOS. This allows us to interpret eV pe ak/2 as the effective hopping parameter t"(U). Our results 
indicate that t"(U) decreases by the renormalization effect due to the interaction. This simple picture 
qualitatively agrees with the SBMFT and the NCA results where the split peaks are located at eV <C t". 
For (U/t, t'/t, t" /t) = (1, 0.6, 0.6) and (2, 0.6, 0.6), not shown in this paper, the splitting cannot be seen 
since the relatively large resonance width V merges the two peaks at ±i"(U) into a single peak. 

The SBMFT and the NCA results are for a fixed energy level in each dot and in the limit U — > oo. If 
we consider the limit of U — >• oo with keeping the particle-hole symmetry, we end up with the state (ii), 
since Tjj <C J c ff and t" <C U/4 always hold. There the pair spin singlet state is frozen and little current 
can flow. Thus, judging from our results we expect that the peak position moves toward V — and, at the 
same time, the current is strongly suppressed by increasing U. 

Now we turn to the left side of the linear conductance peak. G(V) for (t'/t, t" /t) — (0.8, 0.4) shows 
the zero bias peak, Fig|6]right panels. The width of the zero bias peak becomes narrow as increasing U, 
reflecting the formation of the sharp Kondo resonant peak. This behavior is basically the same as in the 
case of the single quantum dot system. In the parameter sets considered here, T is relatively large and 
therefore we see a broad linear response regime. A notable feature in this case is the enhancement of the 
spin correlation function (Si ■ S r ) by the finite bias voltage. The antiferromagnetic spin correlation in the 
steady states turn out to be enhanced compared with the ground state expectation values. Note that the 
states we are focusing on are near the two-Kondo singlet regime (the state (i)) and the spin correlation is 
small since the interdot correlation is effectively suppressed. By forcing the current flow through the double 
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Fig. 6 Current-voltage characteristics J(V), differential conductance G(V) and the expectation values of the spin 
correlation function (Si ■ S r ) for (t'/t,t"/t) = (0.4,0.6) and (0.8,0.4). The TdDMRG calculations have been 
performed for the systems with L = 121. 



dots, the interdot coherence seems to be partly recovered. This does not conform to a naive expectation that 
the finite bias voltage suppresses the electron correlation effects in a similar way as the finite temperature 
does. This enhancement of the interdot correlation by the voltage has similarity to the suppression of the 
double occupancy (ndott^dotj,) by the voltage in the single quantum dot system 1191 1241 . 

4 Conclusions 

In this paper we have reviewed recent developments on the TdDMRG method applied to the nonequilib- 
rium transport through quantum dots. The results for a single dot show clearly that one can obtain reliable 
results on steady currents in systems with correlation effects up to the intermediate coupling regime. 
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In the second part of this paper, we have investigated the nonlinear transport through serially coupled 
double-quantum-dot system. We have accurately simulated the quasi-steady states and discussed how the 
steady currents are affected by the competition between the formation of the two-Kondo singlets (the state 
(i)) and the local singlet states (the pair spin singlet (ii) and the doubly occupied bonding singlet (iii)). 

It has been found that the differential conductance shows splitting of the zero bias peak when the ground 
state is of the local nature, namely near the state (ii) or (iii). The positions of the split peaks in the 
differential conductance are renormalized to smaller values towards as U increases under the condition 
of the particle-hole symmetry. In addition we have shown that in the cases near the state (i) the interdot 
spin correlation can be enhanced by applying a finite voltage. 

The results presented in this paper indicate that it is now possible to perform unbiased numerical cal- 
culations for nonequilibrium steady states of strongly correlated electron systems. The approach based on 
the TdDMRG is from one-dimension and of course a complementary approach is from higher dimensions, 
notably from the infinite dimension. 
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